Depinning transition and thermal fluctuations in the random-field Ising model 
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We analyze the depinning transition of a driven interface in the 3d random-field Ising model (RFIM) 
with quenched disorder by means of Monte Carlo simulations. The interface initially built into the 
system is perpendicular to the [lll]-direction of a simple cubic lattice. We introduce an algorithm 
which is capable of simulating such an interface independent of the considered dimension and time 
scale. This algorithm is applied to the 3d-RFIM to study both the depinning transition and the 
influence of thermal fluctuations on this transition. It turns out that in the RFIM characteristics 
of the depinning transition depend crucially on the existence of overhangs. Our analysis yields 
critical exponents of the interface velocity, the correlation length, and the thermal rounding of the 
transition. We find numerical evidence for a scaling relation for these exponents and the dimension 
d of the system. 
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I. INTRODUCTION 

Driven interfaces in systems with quenched disorder 
display with increasing driving force a transition from a 
phase where no interface motion takes place to a phase 
with a finite interface velocity. This so-called depinning 
transition is caused by a competition of driving force and 
quenched disorder. While the driving force tends to move 
the interface, the motion is hindered by the disorder (see 
e.g. 1). 

Depinning transitions are found in a large variety of 
physical problems, like fluid invasion in porous media [Q, 
depinning of charge density waves or field-driven 
motion of domain walls in ferromagnets |)J . In magnetic 
systems a domain wall separates regions of different spin 
orientations. With the assumption that the correspond- 
ing interface shows properties of an elastic membrane, it 
has been argued || that the depinning of the interface 
can be described by an Edwards- Wilkinson equation || 
with quenched disorder. While the interface motion in a 
system with quenched disorder near the critical threshold 
is theoretically often investigated in the absence of ther- 
mal fluctuations, these fluctuations affect the experimen- 
tal study of the depinning transition The crucial 
point is that energy barriers which are responsible for 
a trapping of the interface in a metastable state at zero 
temperature can always be overcome due to thermal fluc- 
tuations. For driving fields far below the transition field 
this yields a thermally activated creep motion (see M and 
references therein). This behavior changes approaching 
the transition point, where finite temperatures cause a 
rounded depinning transition (for experimental evidence 
see, for instance, Fig. 2 in ||). To describe the depen- 
dence of the interface velocity on driving force and tem- 
perature near the transition point, a scaling ansatz has 
been proposed [Q . This ansatz which is based on an equa- 
tion of motion for sliding charge density waves predicts 
its characteristic velocity to be a power law of temper- 
ature at the critical threshold. This scaling ansatz has 
been shown to be a valid description for the depinning of 
a domain wall in the 2d random-field Ising model (RFIM) 
with quenched disorder [fLOl. 



The outline of our paper is as follows: Sec. || describes 
the RFIM and reflects properties of [lll]-interfaces in 
this model. In Sec. Ill we discuss the depinning tran- 
sition from a microscopic point of view, analyzing the 
mechanisms of interface motion near the depinning tran- 
sition. Also, we determine numerically the exponents 
of the interface velocity and of the correlation length, 
allowing an estimation of the universality class of the 
3e?-RFIM. In Sec. |^ we analyse the influence of tem- 
perature on the depinning transition. By assuming the 
interface velocity to be a generalized homogenous func- 
tion, our analysis is based on applying standard concepts 
of critical equilibrium phenomena. The ansatz allows the 
characterization of the thermal rounding of the depinning 
transition by a critical exponent 5. We determine S for 
the depinning transition in the 3d- RFIM for the first time 
and find numerical evidence for a scaling relation among 
certain critical exponents characterizing this transition. 
This scaling relation also holds in the 2<i-RFIM analyzed 
previously 



II. INTERFACES IN THE RFIM 



We investigate the 3d-RFIM with quenched disorder 
on a simple cubic lattice. The Hamiltonian of the system 
is given by 



H = - J - J2 SiSj - H ~ J2 hi Si ' 



(1) 



(id) 



where the first sum is restricted to nearest-neighbors. 
H denotes the driving field and hi quenched random- 
fields which are uniformly distributed within an interval 
[—A, A]. We carry out Monte Carlo simulations with 
single-spin-flip dynamics and we use transition probabil- 
ities p(Si — > —Si,T), where T denotes the temperature, 
according to a heat-bath- Algorithm (see e. g. |tl| and ref- 
erences therein). At zero temperature these transition 
probabilities reduce to 
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( 1 : sn < 
p(Si ^ -S h 0) = { 1/2 :<5W = (2) 
[ : sn > , 

where = H(—Si) — TL(Si). We investigate three di- 
mensional cubic systems of linear extension from L = 12 
to L = 162. 

An initially flat interface is built into the system sepa- 
rating regions of up- and down spins. The applied field H 
drives the interface. Within the Monte-Carlo simulation 
spins adjacent to the interface flip causing a movement 
of the interface. Also, nucleation may occur, i.e. a spin 
initially parallel to all of its neighbors may turn. Since 
we are interested in the scaling behavior of the inter- 
face motion in the vicinity of the depinning transition, it 
is essential that within the observation time nucleation 
does not occur. The minimum energy needed for iso- 
lated spin flips is 2(zJ — H — A). As long as this quan- 
tity is large as compared to temperature, the time scales 
on which nucleation and interface motion occur are sep- 
arated, and within the observation time no nucleation 
takes place [M. In particular, there is no need to sup- 
press artificially nucleation or isolated spin flips during 
the simulation. 

The analysis of interface motion on simple cubic lat- 
tices considers usually [100]-interfaces. However, investi- 
gating [100]-interfaces in the limit of vanishing disorder 
means that the interface motion is restricted to driving 
fields H/J > z — 2 (see [fL2|). To avoid this, we consider 
[lll]-interfaces which move in the absence of disorder at 
arbitrarily small driving fields increasing the separation 
of time scales for interface motion and nucleation even 
further 0. 

We have found that the most convenient way to imple- 
ment [lll]-interfaces in the numerics is the introduction 




FIG. 1. Periodic images of a moving interface in d = 2. 
Antiperiodic boundary conditions are applied. Black areas 
correspond to Si < and white areas to Si > 0. 
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FIG. 2. Dependence of magnetization reversal processes 
on system parameters in the RFIM at T = 0. The dimension 
enters only through the number z of nearest neighbors. The 
bold line indicates the dependence of the critical field H c on 
the strength of the disorder A. For A < J the critical field 
equals A, while for A > J the bold line is only a sketch. The 
dashed lines indicate the regions where overhangs and island 
growth appear. 

of antiperiodic boundary conditions. This implementa- 
tion is illustrated in Fig. ED. For simplicity, periodic im- 
ages of a snapshot of an interface in d = 2 are shown. As 
can be seen from Fig. |l|, the orientations of up and down 
are exchanged when passing the boundaries of the sys- 
tem. Of course, the exchange of up and down also affects 
the driving field whose sign has to be choosen in an ap- 
propriate manner. Our implementation will work as long 
as the different parts of the interface do not interact. An 
interaction takes place if the interface width, w ~ , is 
of the same magnitude as the typical distance a between 
two neighboring parts of the interface. This distance is 
proportional to the linear extension of the system, a oc L, 
independent of the considered dimension. Our implemen- 
tation is therefore applicable to situations where £ < 1. 
Despite this restriction antiperiodic boundary conditions 
have the advantages that they can be applied to any 
dimension and generalized to other orientations of the 
interface. They are a natural choice for interfaces, be- 
cause the moving interface can be investigated without 
any time limit. This is especially an advantage close to 
the depinning transition, where the critical slowing down 
effect causes large relaxation times O] . 



III. ZERO TEMPERATURE 

In the RFIM with an interface initially built into the 
system, there are in general two magnetization rever- 
sal processes: interface motion and nucleation. Without 
thermal fluctuations the second process does not occur as 
long as (H + A)/ J does not exceed the number z of near- 
est neighbors. The corresponding threshold is shown in 
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FIG. 3. Part of a diagonal interface in d = 2. Different 
orientations of the spins Si are denoted by black circles (fa- 
vored by the driving field) and white circles, respectively. To 
cause a spin flip at the A, B sites, different driving fields are 
necessary (see text). 



Fig. ^ (upper broken curve) . Above this threshold nucle- 
ation processes take place and interfere with the interface 
motion. 

In the following we are interested in the influence of 
overhangs on the value of the critical field H c (A) at which 
the transition takes place. Close to the depinning tran- 
sition, there are two important kinds of spin flips (see 
Fig. g) . All spins of type A with 



S A Sj = will flip if H > H = A , 

(J) 

while the first spin of type B with 



(3) 



Y S B Sj = -2 can flip if H > H_ 2 = 2 J - A . (4) 

(i) 

Here, the sum is taken over nearest-neighbors of A 
and B, respectively. If the strength of disorder A is 
smaller than the exchange energy J, then it follows that 
Hq < H_2- The critical field at which the transition takes 
place is given by H C (A < J) = A . Hence, no overhangs 
occur in the vicinity of the transition point. Taking into 
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FIG. 4. Interface velocity v and its dependence on the 
driving field H for A/ J = 0.7. The depinning transition 
takes place at H — A. The inset shows interface velocities for 
different system sizes L £ {30,42} and ratios A/ J < 1. For 
reasons of clearness not all error-bars are shown. 



account the transition probabilities given by Eq. (|2|) , this 
value of the critical field means that the interface velocity 
depends neither on the driving field nor on the strength 
of disorder as long as no overhangs occur (H < H_2)- 
In particular, in the absence of overhangs the interface 
velocity observed in a disordered system coincides with 
that of a non-disordered system (A = 0). Figure || and 
its inset show numerical data which confirm this scenario 
for the 3<i-RFIM within the error-bars. 

Next we investigate the depinning transition occurring 
in the RFIM for A > J. In this case the transition takes 
place at a certain field H c < A as can be understood 
from the following consideration: For Ho < H-2, not 
all spins of type A can flip if H < A. But because of 
the existence of overhangs, a second growth mechanism 
is possible to the interface: If a spin of type A cannot 
flip due to its large random field hi , an overhang created 
elsewhere can cause an avalanche by which additional 
neighbors of A are flipped. Thus the interface can be kept 
moving. Contrary to the regime A < J, the interface 
motion now is based on the existence of overhangs. Note 
that our considerations do not depend on the dimension 
d of the system, because Eqs. (0) and (Q) are independent 
of d. 

We start examining the regime A > J in the 3d- 
RFIM numerically by investigating the depinning transi- 
tion from below. We analyse the disorder-averaged dis- 
tance (h(t — > oo)) traveled by an initially flat interface 
before pinning occurs. This quantity is closely related to 
the total volume invaded by a growing domain which was 
analyzed in |L4|,[l5). However, while in JT^JT^] the driving 
force is increased step by step to allow for relaxation pro- 
cesses in between, we focus our attention to driving fields 
which remain unchanged during the interface motion. 

Below the depinning transition (h(t — ► oo)) is finite. 
Approaching the transition point with increasing driv- 
ing field, the distance traveled before pinning occurs 
increases and finally diverges at the transition point. 
We assume that in the vicinity of the transition point 
(h(t — > oo)) diverges algebraically, characterized by some 
exponent y, 



(h(t^oo))~(H -H)- 



(5) 



where H c denotes the critical field observed in a system 
of infinite extension. In a finite system with linear dimen- 
sion L finite-size scaling is assumed. The corresponding 
scaling ansatz reads 



(h(t^oo))=LV/ v f\(H-H c )L 1 / v 



(G) 



with f[x) ~ \x\~ v for x — > — oo. Note that (h(t —> oo)) 
also diverges in any finite system which means that f(x) 
should diverge at a finite value of x*. The corresponding 
driving field defines a size dependent critical field H C (L) 
given by (H C (L) - H c ) L x l v = x* . A scaling plot of the 
data according to Eq. @ is shown in Fig. |5|. The diver- 
gence of f(x) occurs at x* « 2.5 showing that in a finite 
system the threshold field is always shifted to fields larger 
that H c . 

The critical exponent of the correlation length parallel 
to the interface is given by \jv = 1.31 ± 0.07 and the 
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FIG. 5. Scaling plot of the distance (/i(i — > oo)) traveled 
by the interface before pinning occurs as a function of the 
driving field H according to Eq. (jfj]). The data collapse yields 
l/v = 1.31 ± 0.07, y/v = 0.98 ± 0.4, and H c = 1.371 ± 0.003. 

critical field turns out to be H c — 1.371 ± 0.03. The 
value of v coincides with ]l5|] , where an [100]-interface in 
the self-affine growth regime corresponding in our case 
to A > J has been investigated. This suggests that the 
behavior of the correlation length at the depinning tran- 
sition does not depend on the orientation of the interface 
in the RFIM. 

In the following we consider the disorder averaged in- 
terface velocity v = (dh/dt) above the transition point 
in the limit of large times. This quantity can be inter- 
preted as the order parameter of the depinning transition. 
Approaching a continuous phase transition the order pa- 
rameter vanishes in leading order according to 



v(H) =A(H - H c 



\0 



(7) 



The corresponding data are shown in Fig. g. The pref- 
actor A is a non-universal constant which can be used to 
compare the results obtained at zero temperature with 
those presented in the next section. Since in the vicinity 
of the depinning transition finite size effects may become 
important, we calculated each interface velocity v(H) in 
systems of different linear extension L. For sufficiently 
large L we observed no significant dependence on the sys- 
tem size from which we concluded that the data shown in 
Fig. |^ correspond within negligible errors to those of the 
limit L — > oo. As can be seen from the data, Eq. ((jj) is 
fulfilled and we obtain A = 0.671 ± 0.03, = 0.66 ± 0.04, 
and H c = 1.37 ±0.01. 

The values of (3 and v obtained from our analysis co- 
incide within the error-bars with those of the Edwards- 
Wilkinson equation with quenched disorder in d = 2 + 1, 
/?ew = 2/3 and vew = 3/4. These values are obtained by 
an e-expansion within a functional renormalization group 
scheme (see (l],[l(|). While the value of /3ew is obtained 
to first order of e, there are arguments that i^ew is exact 
in all orders to e • Taking this into account, our 

results suggest that the depinning transition of a domain 
wall in the 3d-RFIM with quenched disorder is in the 
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FIG. 6. Dependence of the interface velocity v on the 
driving field H in the vicinity of the transition point. Ap- 
proaching the critical field H c , the system size L is increased 
in order to avoid finite-size effects. Fitting the data according 
to Eq. (0) (solid line) yields (3 = 0.66±0.04, H c = 1.37±0.01, 
and A = 0.67 ± 0.03. 



same universality class as the depinning transition of the 
corresponding Edwards- Wilkinson equation. 



IV. FINITE TEMPERATURES 

In this section we study the influence of finite tem- 
peratures on the depinning transition. For T > the 
interface velocity does not vanish for finite driving fields 
since the energy needed to overcome local energy barriers 
is provided by thermal fluctuations at any finite T. This 
results in a rounded depinning transition. The rounding 
can be seen in Fig. ^, where interface velocities for dif- 
ferent driving fields and temperatures are presented. As 
expected, the rounding of the transition increases with 
increasing temperature. Again, we ensured that the in- 
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FIG. 7. Dependence of the interface velocity on the driv- 
ing field for different temperatures, as indicated. The open 
symbols are from Fig. []. The vertical line denotes the critical 
field obtained at T = 0. 
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FIG. 8. Dependence of the interface velocity on H for 
given values of T. The data shown are identical to those in 
Fig. (?) for T > and rescaled according to Eq. (^|). 



terface velocities presented in this and the following Fig- 
ures correspond within negligible errors to those of the 
thermodynamic limit. To analyse the thermal rounding 
of the depinning transition quantitatively, we first note 
that the depinning transition can be described in terms 
of a continuous non-equilibrium phase transition. This 
is suggested by the divergence of the correlation length 
(see determination of v and Fig. ||) and the dependence 
of the interface velocity on the driving field near the tran- 
sition point (Fig. ||). In the standard theory of critical 
phenomena a continuous phase transition is character- 
ized by critical exponents (see for instance |L8| and refer- 
ences therein). Beside j3 describing the field dependence 
of the order parameter and v characterizing the diver- 
gence of the correlation length near the transition point, 
the rounding of a phase transition is characterized by the 
critical exponent 5. In magnetic systems, for instance, (3 
and (5 determine the magnetic equation of state. We now 
apply this approach to the depinning transition by as- 
suming its order parameter to be a generalized homoge- 
nous function of temperature and driving field, 

v [T, H - H c ] = A v [X aT T, \ a » {H - H c )] . (8) 

Choosing A = T~ x l aT we obtain the scaling ansatz 

v(T,H)=T 1 / s f T [(H-H c )T- 1 /P s ], (9) 

with /t(x -»0) = const. In particular, this equation cor- 
responds to the magnetic equation of state ||18|] . From an 
equation of motion of sliding charge density waves a scal- 
ing form corresponding to Eq. (^|) has been obtained [Q . 
Note that contrary to |Q our ansatz which is based on 
Eq. (g) yields no predictions on the values for (3 and S. 

It has been shown previously that Eq. (^) is valid in the 
2d-RFIM with quenched disorder |To| . We have tested 
this scaling ansatz in the present situation for the 3d- 
RFIM with the interface velocities shown in Fig. [?]. As 
can be seen from Fig. ||, the scaling ansatz leads to a 
data collapse for f3 = 0.63 ± 0.07, S = 2.38 ± 0.2, and 
H c = 1.375 ± 0.01. Thus at H = H c the influence of 
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FIG. 9. Dependence of the interface velocity on T for 
given values of H. The data are rescaled according to 
Eq. ([To]). The horizontal line marks the value of A which 
is given by Eq. (^). 



temperature on the interface velocity can be described 
by a power law v ~ T 1 / 5 . To support this value for S we 
can determine 6 from a different scaling function obtained 
from Eq. (§) by choosing A = \H - H^- 1 / 11 ": 



v(T, H) = (H- H c f f H [(H - H c )-^ T] , (10) 

with fu{ x — > 0) = const. This ansatz is valid above the 
transition point and it is closely related to Eq. ([}]). It 
corresponds to a different formulation of the magnetic 
equation of state. Interface velocities rescaled according 
to Eq. © are shown in Fig. ^. One obtains (3 — 0.67 ± 
0.03, S = 2.55 ± 0.37, and H c = 1.37 ± 0.05. This result 
confirms within the error-bars the value of 5 determined 
by Eq. (^). Beside these quantities the data collapse also 
allows a determination of the prefactor A = /h(x — > 0) 
[see Eq. (|)] which turns out to be A = 0.685 ± 0.025. 

The values of A, H c , and (3 found for T > coincide 
within sufficient accuracy with those values obtained at 
T = 0. We have demonstrated that Eqs. (g) and |[(]) 
are valid confirming that the interface velocity is a gen- 
eralized homogenous function in the vicinity of the tran- 
sition point jl9[. Thus, the influence of temperature on 
the depinning transition can be described within well- 
established concepts. 

The knowledge of (3 and 5 allows a test of the scaling 
relation 8 = 2+ 1/(3 proposed by Tang and Stepanow [^0| . 
This scaling relation was shown to be fulfilled in the 2d- 
RFIM @. For (3 sa 0.67 the scaling relation suggests 
S w 3.5 which is not supported by our results. On the 
other hand, standard theory of critical phenomena pre- 
dicts relations among critical exponents. For instance, 
combining the Rushbrooke, the Widom and the hyper- 
scaling relation yields in equilibrium physics 



dv 
J 



1 



fll) 



This scaling relation is valid in dimensions d below the 
upper critical dimension d c due to the restriction of the 
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hyperscaling relation to d < d c . We have tested the 
scaling relation ( |ll| ) with the numerically evaluated ex- 
ponents at the depinning transition and found out that 
both the exponents in the present case d — 3 as well as 
the exponents for d = 2 (i^d ~ 1-0, /?2d ~ 0.33, and 
5 2 d « 5.0; see fuj) fulfill Eq. © within the error-bars. 
Unfortunately however, a firm foundation of this scal- 
ing relation in the present situation for non-equilibrium 
phase transitions is unknown. 

V. CONCLUSION 

We investigated the motion of a driven interface in a 
magnetic system with quenched disorder. To improve 
the efficiency of our numerics we applied antiperiodic 
boundary conditions. These boundary conditions allow 
to investigate the interface motion on any time scale. At 
zero temperature a depinning transition occurs at a finite 
driving field. We discussed the influence of overhangs 
and avalanches on this transition. If the strength of dis- 
order exceeds the coupling constant, the interface motion 
is based on the existence of overhangs. Under these cir- 
cumstances the depinning transition can be characterized 
by critical exponents, both below and above the critical 
threshold. Our results suggest that the depinning transi- 
tion of a domain wall in the 3d-RFIM with quenched dis- 
order and the depinning transition of the corresponding 
Edwards- Wilkinson equation are in the same universality 
class. 

Thermal fluctuations yield a rounded transition. By 
assuming the interface velocity to be a generalized ho- 
mogenous function of temperature and driving field, this 
rounding can be described within a scaling approach. 
The validity of this approach is confirmed by the fact 
that at the threshold field the interface velocity vanishes 
with decreasing temperature according to a power law 
characterized by an exponent S. We have tested a scal- 
ing relation [Eq. (|TT|)1 among different exponents char- 
acterizing the depinning transition and found numerical 
evidence, that the scaling relation is valid both in the 2d- 
and the 3d-RFIM. 
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